Integrated miRNA and transcriptome profiling to explore the molecular determinism of convergent adaptation to corn in two lepidopteran pests of agriculture

BackgroundThe degree to which adaptation to same environment is determined by similar molecular mechanisms, is a topic of broad interest in evolutionary biology, as an indicator of evolutionary predictability. We wished to address if adaptation to the same host plant in phytophagous insects involved related gene expression patterns. We compared sRNA-Seq and RNA-Seq data between two pairs of taxa of Ostrinia and Spodoptera frugiperda sharing maize as host-plant. For the latter, we had previously carried out a reciprocal transplant experiment by feeding of the larvae of the Corn strain (Sf-C) and the Rice strain (Sf-R) on corn versus rice and characterized the mRNA and miRNA responses.ResultsFirst, we predicted the genes encoding miRNA in Ostrinia nubilalis (On) and O. scapulalis (Os). Respectively 67 and 65 known miRNA genes, as well as 196 and 190 novel ones were predicted with Os genome using sncRNAs extracted from whole larvae feeding on corn or mugwort. In On, a read counts analysis showed that 37 (55.22%) known miRNAs and 19 (9.84%) novel miRNAs were differentially expressed (DE) on mugwort compared to corn (in Os, 25 known miRs (38.46%) and 8 novel ones (4.34%)). Between species on corn, 8 (12.5%) known miRNAs and 8 (6.83%) novel ones were DE while only one novel miRNA showed expression variation between species on mugwort. Gene target prediction led to the identification of 2953 unique target genes in On and 2719 in Os, among which 11.6% (344) were DE when comparing species on corn. 1.8% (54) of On miR targets showed expression variation upon a change of host-plant.We found molecular changes matching convergent phenotype, i.e., a set of nine miRNAs that are regulated either according to the host-plant both in On and Sf-C or between them on the same plant, corn. Among DE miR target genes between taxa, 13.7% shared exactly the same annotation between the two pairs of taxa and had function related to insect host-plant interaction.ConclusionThere is some similarity in underlying genetic mechanisms of convergent evolution of two distant Lepidopteran species having adopted corn in their host range, highlighting possible adaptation genes

Molecular evolution of genes in avian genomes

Nam K, Mugal C, Nabholz B, et al. Molecular evolution of genes in avian genomes. Genome Biology. 2010;11(6): R68.Background: Obtaining a draft genome sequence of the zebra finch (Taeniopygia guttata), the second bird genome to be sequenced, provides the necessary resource for whole-genome comparative analysis of gene sequence evolution in a non-mammalian vertebrate lineage. To analyze basic molecular evolutionary processes during avian evolution, and to contrast these with the situation in mammals, we aligned the protein-coding sequences of 8,384 1:1 orthologs of chicken, zebra finch, a lizard and three mammalian species. Results: We found clear differences in the substitution rate at fourfold degenerate sites, being lowest in the ancestral bird lineage, intermediate in the chicken lineage and highest in the zebra finch lineage, possibly reflecting differences in generation time. We identified positively selected and/or rapidly evolving genes in avian lineages and found an overrepresentation of several functional classes, including anion transporter activity, calcium ion binding, cell adhesion and microtubule cytoskeleton. Conclusions: Focusing specifically on genes of neurological interest and genes differentially expressed in the unique vocal control nuclei of the songbird brain, we find a number of positively selected genes, including synaptic receptors. We found no evidence that selection for beneficial alleles is more efficient in regions of high recombination; in fact, there was a weak yet significant negative correlation between ω and recombination rate, which is in the direction predicted by the Hill-Robertson effect if slightly deleterious mutations contribute to protein evolution. These findings set the stage for studies of functional genetics of avian genes

STAT1 and Nmi are downstream targets of Ets-1 transcription factor in MCF-7 human breast cancer cell

AbstractEts-1 is a cellular homologue of the product of the viral ets oncogene of the E26 virus, and it functions as a tissue-specific transcription factor. It plays an important role in cell proliferation, differentiation, lymphoid cell development, transformation, angiogenesis, and apoptosis. Ets-1 controls the expression of critical genes involved in these processes by binding to ets binding sites present in the transcriptional regulatory regions. Here, we transiently overexpressed Ets-1 in MCF-7 and comprehensively searched for potential downstream targets of Ets-1 by cDNA microarray analysis. The expressions of several interferon-related genes including STAT1 and Nmi were augmented by the overexpression of Ets-1. RT-PCR and Western blotting confirmed the increase in the levels of STAT1 and Nmi mRNA and protein. In contrast, Ets-1 siRNA decreased the expression of STAT1 and Nmi proteins. As in our transient transfection experiments, stable overexpression of Ets-1, also increased the protein expression of STAT1 and Nmi in MCF-7 cells. Taken together, our results indicate that STAT1 and Nmi are downstream targets of Ets-1 in MCF-7 human breast cancer cells

Nonlinear Dynamics of Nonsynonymous (dN) and Synonymous (dS) Substitution Rates Affects Inference of Selection

Selection modulates gene sequence evolution in different ways by constraining potential changes of amino acid sequences (purifying selection) or by favoring new and adaptive genetic variants (positive selection). The number of nonsynonymous differences in a pair of protein-coding sequences can be used to quantify the mode and strength of selection. To control for regional variation in substitution rates, the proportionate number of nonsynonymous differences (dN) is divided by the proportionate number of synonymous differences (dS). The resulting ratio (dN/dS) is a widely used indicator for functional divergence to identify particular genes that underwent positive selection. With the ever-growing amount of genome data, summary statistics like mean dN/dS allow gathering information on the mode of evolution for entire species. Both applications hinge on the assumption that dS and mean dS (∼branch length) are neutral and adequately control for variation in substitution rates across genes and across organisms, respectively. We here explore the validity of this assumption using empirical data based on whole-genome protein sequence alignments between human and 15 other vertebrate species and several simulation approaches. We find that dN/dS does not appropriately reflect the action of selection as it is strongly influenced by its denominator (dS). Particularly for closely related taxa, such as human and chimpanzee, dN/dS can be misleading and is not an unadulterated indicator of selection. Instead, we suggest that inconsistencies in the behavior of dN/dS are to be expected and highlight the idea that this behavior may be inherent to taking the ratio of two randomly distributed variables that are nonlinearly correlated. New null hypotheses will be needed to adequately handle these nonlinear dynamics

Genetic structure and insecticide resistance characteristics of fall armyworm populations invading China

The rapid wide‐scale spread of fall armyworm (Spodoptera frugiperda ) has caused serious crop losses globally. However, differences in the genetic background of subpopulations and the mechanisms of rapid adaptation behind the invasion are still not well understood. Here we report the assembly of a 390.38Mb chromosome‐level genome of fall armyworm derived from south‐central Africa using Pacific Bioscience (PacBio) and Hi‐C sequencing technologies, with scaffold N50 of 12.9 Mb and containing 22260 annotated protein‐coding genes. Genome‐wide resequencing of 103 samples and strain identification were conducted to reveal the genetic background of fall armyworm populations in China. Analysis of genes related to pesticide‐ and Bt‐resistance showed that the risk of fall armyworm developing resistance to conventional pesticides is very high. Laboratory bioassay results showed that insects invading China carry resistance to organophosphate and pyrethroid pesticides, but are sensitive to genetically modified maize expressing the Bacillus thuringiensis (Bt) toxin Cry1Ab in field experiments. Additionally, two mitochondrial fragments were found to be inserted into the nuclear genome, with the insertion event occurring after the differentiation of the two strains. This study represents a valuable advance toward improving management strategies for fall armyworm

The genome of a songbird

The zebra finch is an important model organism in several fields with unique relevance to human neuroscience. Like other songbirds, the zebra finch communicates through learned vocalizations, an ability otherwise documented only in humans and a few other animals and lacking in the chickenthe only bird with a sequenced genome until now. Here we present a structural, functional and comparative analysis of the genome sequence of the zebra finch (Taeniopygia guttata), which is a songbird belonging to the large avian order Passeriformes. We find that the overall structures of the genomes are similar in zebra finch and chicken, but they differ in many intrachromosomal rearrangements, lineage-specific gene family expansions, the number of long-terminal-repeat- based retrotransposons, and mechanisms of sex chromosome dosage compensation. We show that song behaviour engages gene regulatory networks in the zebra finch brain, altering the expression of long non-coding RNAs, microRNAs, transcription factors and their targets. We also show evidence for rapid molecular evolution in the songbird lineage of genes that are regulated during song experience. These results indicate an active involvement of the genome in neural processes underlying vocal communication and identify potential genetic substrates for the evolution and regulation of this behaviour. © 2010 Macmillan Publishers Limited. All rights reserved

Recombination Drives Vertebrate Genome Contraction

Selective and/or neutral processes may govern variation in DNA content and, ultimately, genome size. The observation in several organisms of a negative correlation between recombination rate and intron size could be compatible with a neutral model in which recombination is mutagenic for length changes. We used whole-genome data on small insertions and deletions within transposable elements from chicken and zebra finch to demonstrate clear links between recombination rate and a number of attributes of reduced DNA content. Recombination rate was negatively correlated with the length of introns, transposable elements, and intergenic spacer and with the rate of short insertions. Importantly, it was positively correlated with gene density, the rate of short deletions, the deletion bias, and the net change in sequence length. All these observations point at a pattern of more condensed genome structure in regions of high recombination. Based on the observed rates of small insertions and deletions and assuming that these rates are representative for the whole genome, we estimate that the genome of the most recent common ancestor of birds and lizards has lost nearly 20% of its DNA content up until the present. Expansion of transposable elements can counteract the effect of deletions in an equilibrium mutation model; however, since the activity of transposable elements has been low in the avian lineage, the deletion bias is likely to have had a significant effect on genome size evolution in dinosaurs and birds, contributing to the maintenance of a small genome. We also demonstrate that most of the observed correlations between recombination rate and genome contraction parameters are seen in the human genome, including for segregating indel polymorphisms. Our data are compatible with a neutral model in which recombination drives vertebrate genome size evolution and gives no direct support for a role of natural selection in this process

Two Antarctic penguin genomes reveal insights into their evolutionary history and molecular changes related to the Antarctic environment. GigaScience

Abstract Background: Penguins are flightless aquatic birds widely distributed in the Southern Hemisphere. The distinctive morphological and physiological features of penguins allow them to live an aquatic life, and some of them have successfully adapted to the hostile environments in Antarctica. To study the phylogenetic and population history of penguins and the molecular basis of their adaptations to Antarctica, we sequenced the genomes of the two Antarctic dwelling penguin species, the Adélie penguin [Pygoscelis adeliae] and emperor penguin [Aptenodytes forsteri]. Results: Phylogenetic dating suggests that early penguins arose~60 million years ago, coinciding with a period of global warming. Analysis of effective population sizes reveals that the two penguin species experienced population expansions from~1 million years ago to~100 thousand years ago, but responded differently to the climatic cooling of the last glacial period. Comparative genomic analyses with other available avian genomes identified molecular changes in genes related to epidermal structure, phototransduction, lipid metabolism, and forelimb morphology. Conclusions: Our sequencing and initial analyses of the first two penguin genomes provide insights into the timing of penguin origin, fluctuations in effective population sizes of the two penguin species over the past 10 million years, and the potential associations between these biological patterns and global climate change. The molecular changes compared with other avian genomes reflect both shared and diverse adaptations of the two penguin species to the Antarctic environment

Molecular Evolution of the Vertebrate Genome

In this thesis, I studied molecular evolution of the vertebrate genome, focusing on sex chromosomes, protein coding genes, and genome size variation. The evolutionary history of avian sex chromosomes was analyzed by comparing substitution rate among 12 gametologous gene pairs on the Z and W chromosomes. Divergence estimates were distributed into three discrete clusters, evolutionary strata, implying stepwise cessations of recombination. Stratum 3 and stratum 2 are located the intervals 1-11Mb and 16-53Mb on the chicken Z chromosome, respectively. Stratum 1 was located in the middle of stratum 2, suggesting a chromosomal inversion. Using a molecular clock, the estimated times for cessation of recombination between Z and W are 132–150 (stratum 1), 71–99 (stratum 2), and 47–57 (stratum 3) million years ago. Higher divergence rate in the Z chromosome than in autosomes (faster-Z) can be explained by positive selections on recessive alleles in hemizygous females, or by stronger genetic drift due to the smaller effective population size of the Z chromosomes. I found there was no difference in the intensity of the faster-Z effect among male-biased, female-biased, and unbiased genes, as might have been expected under a selection model. This result therefore supports the hypothesis that faster-Z is predominantly due to genetic drift. Next, I analyzed molecular evolution of protein-coding genes in birds. In the comparison of zebra finch, chicken and non-avian outgroups, I found that neutral substitution rate was highest in zebra finch, intermediate in chicken, and lowest in ancestral birds. This difference seems attributable to differences in generation time, ancestral birds being most long-lived. Several functional categories were overrepresented among positively selected genes in avian lineages, such as transporter activity and calcium ion binding. I also found that many genes involved with cognitive processes including vocal learning were positively selected in zebra finch. I also found evidence for Hill-Robertson interference acting against the removal of slightly deleterious mutations at linked loci. Finally, I studied the impact of recombination on genome size variation. I found that highly recombining regions have a more condensed genome structure, including shorter lengths of intron, intergenic spacer, transposable elements and higher gene density. In chicken and zebra finch I found that recombination rate was positively correlated with deletion bias, estimated by sequence comparisons between individual transposable elements (LINEs) and the corresponding master sequences. These observations indicate that the more compact genome structure in highly recombining region is due to a higher rate of sequence loss. Higher deletion bias in autosomes than in sex chromosomes supports this idea. I also found that sequence loss due to the deletion bias can explain nearly 20% of genome size reduction after the split of birds from other reptiles. In human, the recombination rate was positively correlated with the deletion bias estimated from polymorphic indels. These results support the hypothesis that the recombination drives genome contraction via the mutation process